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Superconducting qubit designs vary in complexity from single- and few-junction systems, such as 
the transmon and flux qubits, to the many-junction fluxonium. Here we consider the question of 
wether the many degrees of freedom in the fluxonium circuit can limit the qubit coherence time. 
Such a limitation is in principle possible, due to the interactions between the low-energy, highly 
anharmonic qubit mode and the higher-energy, weakly anharmonic collective modes. We show that 
so long as the coupling of the collective modes with the external electromagnetic environment is suf¬ 
ficiently weaker than the qubit-environment coupling, the qubit dephasing induced by the collective 
modes does not significantly contribute to decoherence. Therefore, the increased complexity of the 
fluxonium qubit does not constitute by itself a major obstacle for its use in quantum computation 
architectures. 

PACS numbers: 74.50.+r, 85.25.Cp 


I. INTRODUCTION 

Starting from the pioneering experiments with Gooper 
pair boxes [T], superconducting qubits have been vastly 
improved [5], thanks in part to better control of the 
qubit environment. For example, 3D transmons [S] ben¬ 
efit from being placed in a cavity that suppress radiative 
losses as well as from their relatively large size that de¬ 
creases the role of surface dielectric losses. Planar trans¬ 
mon variants such as the Xmon [3] have shorter coher¬ 
ence times, but have made possible demonstrations of 
the building blocks of quantum error correction codes 
[ME]- In contrast to the transmon architecture, in which 
a Josephson junction is shunted capacitively, in a flux¬ 
onium qubit the shunt is via an inductance [SIIIQ]. To 
realize in practice a sufficiently large inductance, arrays 
with many Josephson junctions are used. Even in this 
more complicated circuit, long relaxation times have been 
achieved - so long, in fact, to make possible the accurate 
measurement of the phase dependence of the quasiparti¬ 
cle dissipation through a Josephson junction m- Exper¬ 
imentally, the fluxonium coherence time is much shorter 
than the limit imposed by relaxation; in this paper we 
investigate whether the complexity of the circuit, aris¬ 
ing from the use of junction arrays, contributes to this 
shortness. 

The physics of arrays made of identical Josephson junc¬ 
tions has long attracted the interest of theorists and 
experimentalists alike, especially in the context of the 
superconductor-insulator transition controlled by the ra¬ 
tio of charging and Josephson energies [m [T3]. More 
recently, renewed attention has been given to the physics 
of collective modes, whose frequency can be signifi¬ 
cantly lowered below the constituent Josephson junctions 


plasma frequency due to ground capacitances; in particu¬ 
lar, the interplay between collective modes and quantum 
phase slips in rings has been studied pTUlB] . as well as 
the localizing effect of disorder in the junction parame¬ 
ters m Long junction arrays can have a large induc¬ 
tance, proportional to the number of elements, and for 
quantum computation applications such superinductors 
[IB] have been proposed as elements of topologically pro¬ 
tected qubits [TB], for example the O-tt qubit El m- 
In the fluxonium, even in the ideal case of no capaci¬ 
tances to ground, the full symmetry of a ring is broken 
by the presence of a smaller junction. Nonetheless, in 
the ideal case a permutation symmetry is still present; 
the effect of this symmetry and its breaking on the flux¬ 
onium spectrum has been studied in Ref. [22l Aspects of 
the quasiparticle-induced decoherence in the fluxonium 
have been studied theoretically in . Here we focus 

on possible limitation of coherence of the qubit mode due 
to its interaction with the collective modes of the array. 
We will show that both ground capacitances and array- 
junction non-linearities introduce potential decoherence 
channels, but that they are not limiting the fluxonium co¬ 
herence. We nonetheless point out that isolation of the 
collective modes from the electromagnetic environment 
is necessary for long coherence times. 

The paper is organized as follows: in Sec. |ll] we in¬ 
troduce our model for the fluxonium circuit; it includes 
ground and coupling capacitors in addition to the charg¬ 
ing and Josephson energies of each junction. Section [TTl| 
discusses the properties of the odd collective modes, 
which are decoupled from the qubit in the parity and 
time reversal (PT-)symmetric case. The interactions be¬ 
tween even and qubit modes due to ground capacitances 
is the subject of Sec. |IV| In Sec. |V] we consider the ef¬ 
fects of the array-junction non-linearity. In Sec. lYll using 
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the results of the previous two sections we estimate the 
fluxonium dephasing rate due to the interactions with the 
collective modes. The consequences of breaking PT sym¬ 
metry and of placing the coupling capacitors in the array 
are briefly discussed in Sec. Eg and |vmL respectively. 
We summarize our main results in Sec. Ml Numerous 
technical details are presented in Appendices through 
[Hj Throughout the paper we set h= 1. 


II. FLUXONIUM MODEL 


The circuit model for the fluxonium qubit is shown 
in Figg the two ends of an array of W ^ 1 identical 
Josephson junctions (Josephson energy Ej and charging 
energy = e‘^j2Cj) are connected by a so-called phase 
slip junction (Josephson energy Ej and charging energy 
Eq). The loop thus formed is pierced by a magnetic flux 
<i>e- To suppress phase fluctuations in the array, we re¬ 
quire Ej/Eq 1 ; phase slips preferentially take place 
at the phase-slip junction so long as Ej/E^ < Ej/Eq. 
The superconducting islands inside the array have ground 
capacitances (7“ and those at the ends Cg. Identical cou¬ 
pling capacitors of capacitance Cc, used to control the 
system by applying ac voltages +V and —V to them, are 
connected to the islands at the end of the array. The La- 
grangian description of this circuit without the coupling 
capacitors is discussed in detail in Ref.[221 In Appendix [A] 
we briefly summarize (for the paper to be self-contained) 
and generalize (to include the coupling capacitors) the 
relevant parts of that work. 

While there are -|- 1 junctions in the circuit, due to 
charge conservation and flux quantization there are only 
N independent degrees of freedom [22]: the qubit mode 
(j) and the collective modes /r = 1,..., A — 1. In 
terms of these modes, the Lagrangian £5 in the absence 
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FIG. 1. Circuit model for the fluxonium: a small junction 
(top, Josephson energy Ej, charging energy Eq) shunts an 
array (bottom) of many identical junctions. A capacitance 
{Cg or Cg) is present between each superconducting island 
and ground. Additional capacitors Cc, biased at voltages ±U, 
are used to control and read out the qubit, whose spectrum 
depends on the applied flux <l>e. 


of ground and 
= 

Ts = 
Us = 


coupling capacitors takes the form 
Ts-Us, 

— Aj cos {4> -I- ipe) 

~^J ^ cos — + ^ , 

m \_ fi 


( 1 ) 

( 2 ) 

(3) 


where the phase bias (pe is due to the externally applied 
magnetic flux, ipe = 27r$e/‘ho with $0 the flux quantum, 
and the qubit mode charging energy Eq is 


1 _ 1 1 


(4) 


note that due to the last term in the above equation, Eq 
is smaller than the phase slip junction charging energy 
Eq. Hereinafter, sums over index m run from 1 to A and 
those over Greek indices such as p from 1 to A — 1. To 
write the Lagrangian in the given form, the matrix W^m 
must satisfy J2m Wf,mW„m = and Y.m = 0 ; 

for concrete calculations, we will use for the form 

suggested in Ref. |22| 


TiA 7r/r(m-l/2) 

W^m - Y ^ cos ^ -. (5) 

It can be shown [22] that £5 is symmetric under the 
action of the symmetric group Sat. 

Due to the large ratio Ej/Eq, fluctuations of are 
small; for fluctuations in (j) small compared to ttA, we 
can then expand the last term in Us to quadratic order 
to find (up to a constant term) 

Us^U = -E’} cos (<^ + Pe) + i A} ^ ( 6 ) 

with El = Ej/N. At this lowest order, the qubit mode 
(j) and collective modes do not interact (so long as 
we neglect ground and coupling capacitors), and the La¬ 
grangian 


Cu=Ts-U (7) 

is symmetric under the unitary group U(A — 1). We note 
that the symmetry under permutations ensures [ 22 ] that 
the anharmonic terms that we neglect cannot couple any 
state given by the direct product by a qubit eigenstate 
and a collective modes singly excited state to a state 
which is the direct product between any qubit state and 
the collective modes ground state. One of the main ob¬ 
jectives of this work is to understand the interactions in¬ 
duced by the presence of ground and coupling capacitors; 
we will show, for example, that the coupling between the 
states just discussed is in general present, albeit weak. 
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Ground and coupling capacitors modify the kinetic en¬ 
ergy part of the Lagrangian by the addition of the term 
To given by 


Tg = 


1 


16 


+ 2 ^ -I- ^ 

^ ITU 

where the symmetric matrix G has entries 




Goo — 
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Go^ = - 


2N + 
3 N 
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-A 


2E- 


1 






( 8 ) 

(9) 

( 10 ) 

( 11 ) 


Here the energy scales are Eg = l2Gg and Et = 
l2Gt, where the total capacitance Gt is 

Gt = 2 {Cl + C,) +{N- 1)C“ . (12) 

The dimensionless parameter 0 < A < 1 is defined as 

(iV-l)C“ 


A = 


Gt 


and we introduced the short-hand notations 

. TT/i 

Su = sm , 

M 2N 

TT/r 

c„ = cos , 

M 2N 


i-{-iy 


(13) 

(14) 

(15) 

(16) 

The coupling capacitors enable control of the circuit 
via external ac voltages; the coupling Lagrangian Cy 
takes the form 

= (17) 

where we assumed that the voltages applied to the ca¬ 
pacitors are equal in magnitude V and opposite in sign 
and E'^ = e^/2Cc. 

From the above definitions, it is evident that the collec¬ 
tive modes with even index ^ interact only with the qubit 
mode, whereas the odd modes interact only among them¬ 
selves. Therefore the total approximate [55] Lagrangian 
£ separates into a sum of even and odd sectors: 


III. COLLECTIVE MODES: ODD SECTOR 

The collective modes odd under PT-symmetry are gov¬ 
erned by the Lagrangian 


^o = To--E:}j2y 

p=i 


No 


(19) 


where to simplify the notation we introduce = ^ 2 p-i 
and the number of odd modes No = [N/2\ equal to the 
integer part of N/2. In the kinetic energy term, we sep¬ 
arate a purely diagonal term Td independent of A , and 
a term Tx proportional to A: 


To — Td + XTx , 
N, 


1 


^‘'=16^ 


Tx = --. 


p=i 

1 
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E: 
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No 
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■S 2 p-l 


Sp 


C2<T-lC2p-l 


If the number of junctions is not too large, 


( 20 ) 

( 21 ) 

aCp- ( 22 ) 


(23) 


the effect of the ground capacitances can be treated per- 
turbatively for all modes and all values of A; if the ar¬ 
ray is short compared to the “screening length” Ng, the 
ground capacitances hardly affect the energy spectrum 
of the modes m- The condition ( |^ is not satisfied 
in current experiments, see Table |lj However, as the or¬ 
der number a of the modes increases, the effect rapidly 
diminishes m-, moreover, the off-diagonal part is pro¬ 
portional to A, which is typically somewhat smaller than 
unity. This points to the more general viability of a per¬ 
turbative approach than suggested by the condition (23) 
above. Formally, diagonalization of the kinetic energy 
matrix can be obtained by a rotation that eliminates the 
off-diagonal terms [see Appendix [C| note that since the 
potential energy term in Eq. (19) is quadratic and pro¬ 


portional to the identity matrix, any rotation leaves it 
unchanged]; up to second order in A we find: 


£ = £[/ +Tg + Cv = Le + Co ■ (18) 

This separation is valid as long as parity P and time re¬ 
versal T symmetries are preserved |22| : indeed, the qubit 
mode and the collective modes with even index ^ are 
even under PT symmetry, while modes with odd index 
are odd. We restrict our attention to the PT-symmetric 
case for most of the paper, but we discuss the conse¬ 
quences of breaking this symmetry in Sec. |VH| In the 
next section we focus on the odd sector. Throughout 
the paper we will present examples calculated with the 
two parameter set given in Table |l| as explained in Ap¬ 
pendix]^ the parameters are chosen as to reflect realistic 
experimental values UM- 


No 




p—1 '-'^P 


(24) 



N 

E^j 

Eg 


Eh 


El 

Eh 

A 

Ns 

set 1 

43 

26.0 

1.24 

8.93 

3.60 

194 

6 

24.2 

0.34 

39 

set 2 

95 

48.3 

1.01 

10.2 

4.78 

484 

807 

12.1 

0.54 

69 


TABLE 1. Fluxonium parameters used for numerical calcu¬ 
lations throughout the paper (see also Appendix 0. Ener¬ 
gies are given in GHz, while the number of jun ctions N, A 
[Eq. ( |13[ )], and the screening length Ns [Eq. ( |23[ )] are dimen¬ 
sionless. 
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where the effective charging energy for the new, rotated 
modes Qp is 


^C,p 
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4B; 
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‘-2p-l 


+ A 


g L^2p-1 
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N{N-1) 


E' 


^(^- 1 ) 4-1 

„2 „2 

2p-ll-2o--l 


t/p 2 p 


-1^2cr-l 


'S2p-l'®2cr-l J 


(25) 


Comparison with numerical diagonalization for a range 
of experimentally relevant parameters shows that this 
formula is accurate to a few percent even for the lower 
energy modes (and much better accuracy for the high- 
energy modes), and that the term proportional to A^ does 
not contribute significantly to the effective charging en¬ 
ergy, thus validating our perturbative approach. With 
the Lagrangian in diagonal form, it is straightforward to 
obtain the energy spectrum of the modes: 

4 = ■ ( 26 ) 

In Sec. I VII A| we will briefly compare the calculated spec¬ 
trum with recent experiments [181128] . In the next section 
we turn our attention to the even sector. 


IV. EVEN SECTOR AND QUBIT MODE 


The qubit mode 4> belongs to the even sector, and in 
the presence of capacitance to ground in the array it in¬ 
teracts with all the even modes: 


We 


Ce — + Te - + Cy 


1 1 . 


p=l 

We 


(27) 




E C2p 


9 p=l 2p 


Here we define rjp = ^ 2 p and the number of even modes is 
iVe = [{N — 1)/2J. The kinetic energy of the even modes 
has a simple diagonal form 


1 


p=i 


1 1 1 
^^9 ^2p 


iVe 


^4 16 E 

p—1 ^■,p 


and the qubit Lagrangian is 
1 


Erh — 


16E. 


c 


+ EjCOs{(l) + ipe) - -El(, 


(29) 


where the effective qubit charging energy is given by 


—1 f 2 N + 1 \ 
E^~ E^^ 4,Et\ 3 N ) ' 


(30) 


The Lagrangian has a simple structure, describ¬ 
ing a set of independent harmonic oscillators interacting 
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FIG. 2. (Color online) Qubit-collective modes coupling con¬ 
stants Qp, Eq. (351, calculated using parameter set 1 in Ta- 
ble|T] (filled circles) and 2 (empty cricles). Different horizontal 
scales are used for the two set, as indicated by the arrows, and 
hence gp for set 2 is larger than the corresponding coupling 
constant for set 1. Note that in absolute value the coupling 
strength between qubit and the few lowest even modes is as 
strong as or even stronger than the typical qubit-cavity cou¬ 
pling [sno]. 


with the qubit mode. To this Lagrangian, however, corre¬ 
sponds an Hamiltonian Hg in which all degrees of freedom 
interact among themselves [see Appendix]^ due the non 
diagonal form of the kinetic energy - cf. the last term in 
Eq. ( |2^ . The condition in Eq. ( |M| ) is again sufficient to 
enable perturbative treatment of these interactions but, 
as mentioned above, it is not experimentally satisfied. In 
contrast, the generally weaker condition 

(31) 

is satisfied in current experiments, with the right hand 
side being about 6400 (10600) for parameter set 1 (2). 
This conditions enables us to make substantial simpli¬ 
fications, with the approximate Hamiltonian taking the 
form 

We 

H, = H^E^HpE LTint + Hy , (32) 

p=i 

= AE^pI - 4 cos + Pe) + \el4>‘^ , (33) 

Hp = AE^CpPI + \eWp , (34) 

cc ^ 4 Et,E^c,pC2p ,,,, 

Hint — 2_^ 9pPpP4> > 9p — E°- 

E't’ 1 

Hv =-‘2^P<j>eV-'^—gpPpeV. (36) 

c p=i ^ c 

Here Pp is the momentum conjugate to rjp. The above 
expression for the Hamiltonian is one of the main re¬ 
sults of this paper: it contains the leading interaction 
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terms between the qubit mode and the collective modes 
of the junctions forming the superinductance due to ca¬ 
pacitance to ground in the array. The coupling constant 
gp is proportional to the array capacitance to ground and 
monotonically decreases with p, see Fig. therefore the 
higher collective modes couple more weakly to the qubit. 
Similarly, since gp determines also the coupling of the 
collective modes to the ac voltage V^ the higher modes 
are more weakly coupled to V than the lower ones; more¬ 
over, the low-energy ones are more weakly coupled to V 
than the qubit mode. As we will see below, this implies 
a lower decay rate for the collective modes compared to 
the qubit. 



A. Dispersive shifts 


To further study the effect on the qubit of the interac¬ 
tion with the collective modes, we perform a Schrieffer- 
Wolff transformation and project the Hamiltonian Hg, 
Eq. (32), onto the qubit subspace to find the effective 


Hamiltonian 


HeB = 


wio(/) 


iVe 


+'^[^P + Xp{f)^lalap, (37) 

p=i 


where is a Pauli matrix in the qubit subspace. 


^p ~ 


(38) 


is the harmonic oscillator frequency of the even collec¬ 
tive mode p, and a^, Op the creation and annihilation 
operators. As indicated by the presence of the parame¬ 
ter / = (pe/27r, the qubit frequency wio and the ac Stark 
shifts Xp depend on flux through the loop, and we ne¬ 
glect for the moment the coupling to external bias given 
by Hy- Here with the frequency wio we indicate the en¬ 
ergy difference between the two lowest eigenstates of the 
Hamiltonian Eq. (33); below we will discuss a small 
renormalization Jwio mAhe qubit frequency due to the 
interaction with the collective modes. 

The dispersive shifts Xp depend on matrix elements 
of charge operator that involve all the eigenstates |Z) 
with energy e; {I = 0,1, 2,...) of Hamiltonian |29) : 


1 


Xp = 


E' 




gl l(ob<>|i)r 


2uj 


10 


W?n - 


l>2 

l>2 


^10 


wFo - ^p 




- ^p 


(39) 


with uJij = ei — ej. While in the case of the transmon [30] 
selection rules and low anharmonicity enable the analyt¬ 
ical calculation of the dispersive shift, for the fluxonium 
this is not possible, so we resort to numerical estimates. 


FIG. 3. (Color online) Main panel: matrix elements squared 
|(0|p,^|Z)P (filled circles) and |(l|p,/,|Z)p (empty circles) at / = 
0.35 for Z < 20 calculated using parameter set 1 [see Table [I|. 
Inset: same as main panel but for parameter set 2. We stress 
that the overall decrease of the matrix elements with Z is valid 
for any flux /, not just for the particular value used here as 
an example. 


Fast convergence of the sums in the above equation is 
ensured by the rapid decrease of the matrix elements 
(O|p 0 |Z) and (l|p 0 |Z) as I increases, see Fig.further aid¬ 
ing the convergence is the approximately linear increase 
of the energy of the states with slope ~ (SE^Al)^/^. The 
reason for the decay of the matrix elements is the follow¬ 
ing: the low-lying states are localized near (j)g, while the 
high-energy states are to a good approximation the eigen¬ 
states of the harmonic oscillator obtained by neglecting 
the Josephson term in (this also explain the linear 
increase in their energy). Since these high-energy states 
display oscillations of small magnitude at the center of 
the potential, the overlap of their derivative with the low- 
lying states is small, and the increase of the number of 
oscillations with I causes the decrease of the matrix ele¬ 
ments. 


The dispersive shifts calculated via Eq. (39) diverge as 
the energy differences ujiq or wn approach one of the col¬ 
lective modes frequencies uip. This divergence, however, 
only signals the breakdown of the perturbative calcula¬ 
tion near such resonant conditions: the actual dispersive 
shift is limited in magnitude by the coupling constant 
gp, so validity of the perturbative approach is given by 
the condition \xp\ AC gp |3T|. Despite this limitation, 
Eq. (39) correctly estimates the dispersive shifts at most 
flux values. We show in Figs. [^ and the flux depen¬ 
dence of dispersive shifts xi X 2 ^ respectively. We see 
that away from resonances the shifts for set 2 are usually 
larger, as expected due to the larger coupling strengths 
[cf. Fig. i- Moreover, due to the higher energy of the 
collective mode involved, the shifts X 2 displays a richer 
resonance structure. We will comment on the effect of 
these shifts on qubit coherence in Sec. |VI| Next, we con¬ 
sider the effect of the qubit-mode coupling on the qubit 
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FIG. 4. (Color online) Dispersive shift xi [Eq. (39l] as func¬ 
tion of external flux for parameter set 1 (thin line) and 2 
(thick line). 
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FIG. 5. (Color online) Dispersive shift X 2 [Eq. @] as func¬ 
tion of external flux for parameter set 1 (thin line) and 2 
(thick line). 


frequency in the absence of excitations of the modes. 


FIG. 6. (Color online) Contributions 5a)io,i (thick line) and 
SijJio ,2 (thin, dashed) [Eq. (41l] to the qubit frequency renor¬ 
malization dcuio for parameter set 1 (inset) and 2 (main 
panel). For comparison, we note that the unrenormalized fre¬ 
quency ui\o approximately varies between 9 GHz (8.2 GHz) 
at zero flux and 0.33 GHz (0.64 GHz) at half flux quantum 
for parameter set 1 (set 2). 


with 


X M 2 


2 V8^C,P " 


|( 0 b^|I)|^^^ 


+EKob^ior 

l>2 


1 


ai/0 + 
1 


(41) 


l>2 


^11 + ti-'p 


there is one im- 
wio, there are no 


While this formula resemble Eq. 
portant difference: so long as ujp > 
divergences in Eq. (41). In Fig. we plot the first two 
largest contributions to dwio, namely Jwio.i and i5a;io,2- 
The former is generally much larger (in absolute value) 
than the latter, due to the stronger coupling [cf. Fig.[^. 
Note that even at half flux quantum, where wio has a 
minimum [Mill], the correction is at most a few per¬ 
cent of WlQ. 


C. Purcell rate 


B. Qubit frequency renormalization 


In the preceding section we studied the change in the 
qubit frequency when a collective mode is excited, but 
the interaction term Hmt of Eq. (35) modifies the qubit 


frequency even in the absence of collective mode excita¬ 
tions, Wio —t Wio -fJwio. This frequency correction dwio, 
arising from Lamb-type energy level shifts, is given by 


dwio = Sui 


10,p 


(40) 


So far we have considered the system to be capacitively 
coupled to external voltage sources. In practical realiza¬ 
tions of circuit QED experiments, this coupling is to a 
mode of a cavity; this can be accounted for by replacing 


V V(c^ + c) 


(42) 


in the coupling Hamiltonian Hy-, Eq. (36). Here parame¬ 


ter V accounts for the strength of the electric field at the 
qubit position as well as for the geometry of the cavity- 
qubit system, while (c) are creation (annihilation) op¬ 
erators for photons in the cavity. As it is customary, to 
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include the cavity and its coupling to an external bath 
of harmonic oscillators, we add to in Eq. (32) the 
following Hamiltonian H^h 


Hcb = uJcC^c + uJkblbk 


Xk 


, (43) 


where w. 


is the cavity frequency, b], (bk) are creation 


(annihilation) operators for bath excitations with energy 
uik, and Xk the coupling strengths between cavity and 
bath modes. 

Within this model, one can calculate using Fermi’s 
golden ruled the so-called Purcell rate Hq for the qubit; 
that is, the decay rate of the qubit excited state by emis¬ 
sion of a photon into the bath (mediated by the cavity): 


«?(/) = K 


glU) 


i^ioif) - ^cY 


(44) 


where n is the inverse lifetime of a photon in the cavity, 
as determined by the cavity-bath couplings Xk |30j . and 


9 <! = 2eH^(0|p,^|l). 


(45) 


This coupling constant depends on flux via the qubit 
states. The above expression for Kq is valid in the dis¬ 
persive regime gq <C jwio — Wc|. The similar calculation 
for the collective modes gives their decay rate as 


1 


Ko = K- 


{Ojp - UJc) 




eP 


1 




(46) 


where ip = {8EQpjEjY^‘^ and again we have assumed 
that the factor multiplying k is small compared to unity. 
Note that the term in square brackets is a decreasing 
function of p; therefore, so long as all the modes have 
frequency above the cavity one (wi > ujY), the decay rate 
of the collective modes decreases with p. If the qubit 
lifetime is limited by Purcell relaxation, we can then es¬ 
timate the modes’ lifetimes by elim inating the unknown 
quantity V from Eqs. (44) and (46). In fact, the ratio ra¬ 


tio Kp/Kq{f) is also independent of k and is determined 
by the circuit properties and the cavity frequency; we plot 
Kp/Kq{0) in Fig. 0 for the parameters given in Table 1^ 
with = 8.18 GHz for set 1 and = 8.89 GHz 

for set 2 [see Refs.and [TTl respectively]. 


o 
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FIG. 7. Ratio KpjKqiQ) between collective mode and qubit 
Purcell decay rates [Eqs. ( |46[ ) and ( |44[ ), respectively] for the 
ten lowest even modes. Filled circles: parameter set 1 in 
Table 12 empty circles: set 2. 


begin with, we split the Josephson energy of each junc¬ 
tion into two contributions using the identity 


cos 




COS — cos 
N 


— sin — sin 
N 






(47) 


The sine product term, as we discuss below, generates 
two qubit-collective mode interaction terms that we will 
denote with and . The product of the two 
cosines can be rewritten as: 


cos — cos 
N 


cos 


E 
E 




-c, 


Om 


(48) 


+ I cos — - 1 1 I cos 


E 


L M 


-C, 


Om I ) 


where 


V. NON-LINEARITY OF THE ARRAY 
JUNCTIONS 


i'Om — ( COS 


E 


(49) 


By expanding the last term in Eq. (§ to quadratic or¬ 
der to obtain the approximate potential energy term in 
Eq. (§, we have treated the Josephson junctions in the 
array as linear elements. However, the cosine in Eq. @ 
includes their non-linear properties, and in this section 
we account perturbatively for these non-linearities. To 


is the expectation value of the operator inside the angu¬ 
lar brackets in the ground state of the collective modes. 
As detailed in the next section, the first term in the right 
hand side of Eq. (481 gives rise to the qubit mode effec¬ 
tive potential while the last term is a qubit-collective 
modes interaction contribution, denoted with ■ The 



























third term is a constant that can be neglected. The sec¬ 
ond term gives, to lowest order, the harmonic potential 
energy of the collective modes. 


E 


cos 


E 




(50) 


The higher order terms in the expansion neglected here 
leads to interactions among the collective modes that do 
not affect the qubit directly. We do not consider such 
interactions from now on, and write the potential energy 
Us in the approximate form 








<Pi 


(51) 


1=1 


U) 


with the potentials and 


specified in what follows. 


A. Qubit effective potential 


Keeping only the first term in the right hand side of 
Eq. (48), from Eq. ^ we find 


= -Ej cos {(j) -I- ipe) - El 


cos — 
N 


(52) 


with Co = J2m^0m/N. Upon expansion of the term 
in square bracket (valid for \(j)\ <C ttN) and assuming 
Cq « 1, we recover the quadratic inductive energy term 
ELti^/2 of Eq. (^. The full qubit potential, however, 
contains small additional non-linearities originating from 
the higher-order terms of the expansion. Here we do not 
consider these terms further, as they are suppressed by 
factors of the form 1/N^^ (j = 1, 2, 3,...), but we show 
that in general Co < 1; therefore, the actual inductive 
energy 


El = ElCq 


(53) 


is smaller than what the simple expression El = Ej/N 
suggests. 

The expectation value entering Comi Eq- can be 

readily obtained from the known matrix elements for 
the harmonic oscillator, see for example Appendix D of 
Ref. [2S1 Here we have to remember that in the odd sec¬ 
tor a rotation from the original modes Cp to independent 
modes (p is necessary, see Sec. HI this is acc ompl ished 
via the orthogonal matrix Ap^, denned in Eq. (C5). We 
thus arrive at 




N 


-jYl Wp^W^^ArpArX 




, (54) 


where 


ir = i8Ec,rlE'})^/^ 


(55) 


is the oscillator length for mode r, with Ec,t given in 
Eqs. and ([2^ for odd and even modes, respec¬ 


tively. Clearly Cq < 1 so long as at least one oscillator 
length is finite. We can also find a lower bound (and 
rough estimate) for Cq by noting that ir < £o, where 
£o = {8 Eq/E j)^^‘^ is the oscillator length in the absence 
of capacitance to ground in the array; note that we typ¬ 
ically have ^0^1; cf. Table |l] Then, using the identities 
Er = 1 and = {N- 1)/N, we find: 


Co > exp 


-r 

tn 


,N-l 

47V 


(56) 


The expansion to lowest order in £q of this formula agrees 
with the expression for the reduction of El reported in 
Ref. [22l Our result shows that that expression generally 
overestimates the suppression of the inductive energy. 


(2) 

B. Quadratic interaction 


We now consider the leading order contribution to 
the potential energy originating from the last term in 
Eq. p8| . By expanding the term dependent on the 
collective mode coordinates and introducing the cre- 
ation/annhilation operators via ^p = £p{ap -I- ap/-\/2 we 
arrive at 


^ (1 - cos ^ ) E W (^OpOp + a)pa\ + UpUp) . 


N 

(57) 

Note that in the absence of the array ground capacitances 
(so that £p = £o) this interaction term is invariant un¬ 
der orthogonal transformations belonging to the group 
0{N — 1) - the U{N — 1) symmetry of approximate La- 
grangian Cu in Eq. Q is only partially broken (the first 
term in brackets actually fully preserves U(A^ — 1) sym¬ 
metry) . 

In Eq. (571 we can distinguish two contributions. First, 
there are terms which are proportional to each collective 
modes number operator Up = aj^Up; to lowest order in 
1/N these terms are 


Us El 



2N 


^ ^ ^p®p®p ) 


(58) 


and they give a dependence of the inductive energy El 
on the occupation of the collective modes. Since changes 
in El lead to variations of the qubit frequency, this de¬ 
pendence can be interpreted as a dispersive shift Xp^'" ■ 


SEl ^^10 ^ EE/2 

2 V QEl y 2 a ^ ■ 


(59) 


Over a broad range of fluxes, except near half-integer 
multiples of the flux quantum, the qubit frequency is ap¬ 
proximately proportional to the inductive energy |23| . 


Wio 


{2EfEL 


1 

2 


(60) 






















9 



FIG. 8. Normalized derivative diox^jOE l calculated for the 
parameters in Table [T] Solid line: set 1; dashed line: set 2. 
Note the minima a,t f — 0 and / = 1/2. 


SO that El (duJio/dEL) « wiq. This approximate rela¬ 
tion translates at all fluxes in an upper bound for the 
dispersive shifts: 


sEl <r ^10 a 
-- 4N ^ ' 


(61) 


This bound shows that the dispersive shifts lead to rela¬ 
tive changes of order 1/fV, in the qubit frequency. Inter¬ 
estingly, the derivative duJio/dEL and hence the disper¬ 
sive shifts have minima at half-integer multiples of the 
flux quantum, see Fig. Therefore the dephasing in¬ 
duced by Us El is suppressed at these “sweet spots”, sim¬ 
ilar to the suppression of dephasing by flux noise; this is 
not surprising, since at leading order the flux and induc¬ 
tive energy affect the qubit frequency in the same way, 
see Eq. ( |60[ ). However, for typical experimental parame¬ 
ters the reduction is less than one order of magnitude. 


( 2 ) 


comes from 
the effect 
Schrieffer- 


IV 


The second type of contribution in 
terms of the form -I- fli . As in Sec. 
of such terms can be studied by performing a 
Wolff transformation, as detailed in Appendix]^ Here we 
simply note that, since they involve the virtual exchange 
of two collective mode excitations rather than one, the 
resulting dispersive shifts are generally smaller than 
of Eq. (591 and can therefore be neglected. 


C. Linear interaction U 


( 1 ) 


In this subsection and the next one, we focus on the 


perturbative treatment of the last term in Eq. (471, ob¬ 


tained by expanding the sine with argument the collective 
modes coordinates. The contribution to the potential en¬ 
ergy from the linear term in this expansion vanishes by 
construction, due to the property = 0. As we 

show in Appendix the third order term gives rise to 
two types of interactions, one of them being a linear in¬ 


teraction term of the form 


[/(i) _ 




Y^gp i^NSin - 
1 


(62) 


3p - 2(2iV)3/2-^“ 


with rjp and defined in the text after Eq. (27). Since it 
couples the qubit mode with even collective modes only, 
this interaction preserves PT-symmetry. Moreover, in 
the Sat symmetric case (i.e., neglecting ground capaci¬ 
tances, so that for all g) this term is absent, in 

agreement with the group-theoretical analysis of Ref. 
The coupling constant Qp decreases with the collective 
mode index p, albeit more slowly than pp in Eq. (35) for 
low p, while for large index p < Nf, we find 


<E 

9p 


32NE^ 


(63) 


where a;“ = ^SEjEq is the array junction plasma fre¬ 
quency. A (loose) upper bound for gi is given by 


9i < 


S'? 


^0 


2 (2V)3/2 2 (2iV)3/2 ■ 


(64) 


As done for the interaction term i7int in Eq. (351, the 


effect of on the qubit can be more easily studied by 
performing a Schrieffer-Wolff transformation leading to 
the additional dispersive shift X^p\f)' 


( 1 ) 1 ^^C,p ~2 


+E 

l>2 

-E 

l>2 


iVsin^lIl) 


iVsin^llO 


2uj 


10 


- ^p 




(1| (iVsin^ ) 10 


Wil 




(65) 

A few comments are in order: first, since the low-lying 
states are localized (so that relevant values of 0 are at 
most of order 27r) and N is large, a good approximation 
for the matrix elements is obtained with the substitu¬ 
tion N sin (j)/N —> (f). Second, within this approximation, 
numerical calculation of x^p'^ does not require much ad¬ 
ditional computation compared to that of Xp, Eq. (39), 
thanks to the identity [35] 




= (SE^^ \{m\p^\l)f . (66) 


Third, the pole structure in Eq. (65) is the same as that 


of Xp and the matrix elements again becomes smaller as 
I increases, see Fig. and Eq. (66). Finally, for typical 
experimental parameters, the coupling constants Pp are 
at least two orders of magnitude smaller than pp [see 
Fig. 10 , so we expect the dispersive shifts Xp^^ to have 


negligible effect - we will return to this point in Sec. VI 
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FIG. 9. (Color online) Main panel: matrix elements squared 
|{0|(()|Z)P (filled circles) and |(1|<()|Z)P (empty circles) at / = 
0.35 for Z < 20 calculated using parameter set 1 [see Table |I]. 
Inset: same as main panel but for parameter set 2. 



modes ii ^ = v and three modes \i ^ ^ v. Note that the 
index structure ensures that at least one index is even 
and the remaining two indices have the same parity; this 
shows that PT-symmetry is preserved. 

As a consequence of the presence of three creation- 
annihilation operators in Eq. (67), a description in terms 
of an effective Hamiltonian would involve terms with 
products of up to three number operators (see also the 
discussion of the interaction in Appendixj^. Rather 
than attempting such a complicated description here, we 
consider the case in which the occupation probability of 
each mode is sufficiently small that we can neglect the 
possibility of having two or more excitations in a mode 
or two or more modes being excited at the same time; 
this requires the occupation probability to be small com¬ 
pared to 1/Af, see Appendix [f| In other words, we only 
consider the possibility that no more than one collective 
mode is excited at any given time. We can then calcu¬ 
late the change Awio^^ in qubit frequency from when the 
collective modes are in their ground state |0) (a^|0) = 0 
for any /r) to when one of the collective modes is excited, 
i.e., in state |1^) = oj) |0). Such a frequency change re¬ 
sembles the dispersive shifts discussed so far, although 
those are valid for multiple excitations in each mode. 
The perturbative calculation of the frequency change is 
detailed in Appendix there we also show that the fre¬ 
quency change is smaller than the dispersive shift 
in Eq. ( [^ and can therefore be neglected. In the next 
section we explore some effects of the collective modes 
dispersive shifts on the qubit. 


VI. QUBIT DEPHASING 


FIG. 10. Goupling constants ratio gpjgp for parameter set 
1 (inset, filled circles) and 2 (main panel, empty circles) for 
all the even modes. Horizontal dashed lines are given by the 
right hand side of Eq. (631. 


D. Multi-mode interaction U' 


(3) 




All the interactions discussed so far involve the qubit 
mode and a single collective mode. This is not the case 

(3) 

for the interaction term , which involves up to three 
collective modes: 




=-^ sm — 

24\/iV N 




3 ^ 

N-1 

p-\-v>N 

X ■^27V-ri-i/H.C.^ 

(67) 

where H.c. denotes the Hermitian conjugate. The first 
term, for example, contains creation operators of two 


As it is well known, the last term in the effective 


Hamiltonian Eq. (371 can be interpreted as a shift in the 


qubit frequency dependent on the states of the collective 
modes; more generally, we write 


wio =^io + 2y^x 




( 68 ) 


where rip^ is the occupation number of mode fi and we 
use X/i to denote the total dispersive shift of that mode. 
For even /r = 2p it is given by 


X2p — Xp + Xp^^ + 


(69) 


while Xfj. = Xp ^ for /i odd. Therefore for even modes 
the total shift includes the contributions from Eq. (39) 
and (65) together with that in Eq. (59). However, we 
note that for almost all values of flux [35] we have 
{gp/gpY and the latter quantity is < 10“® for 


IXp^VxpI 


experimentally relevant parameters, see Fig. |10[ there¬ 
fore, we neglect Xp^^ from now on. As for relative impor¬ 


tance of the contributions Xp and X 2 ^^j note that 
while the former quickly decreases in magnitude with 
increasing index p, the latter slowly increases with p. 
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Therefore we can expect that while it may be necessary 
to keep Xp for the low index modes, could be the 

only relevant contribution for higher index modes. To 
identify which modes have low index in this sense, we re¬ 
mind that in the dispersive regime \xp\ ^ ffp, so we can 


certainly neglect Xp if 9p < IX 2 p^l- Unfortunately this 
latter condition is satisfied at any flux only for high index 
modes, so in general we must keep both Xp and X 2 ^^ fo'^ 


quantitative estimates - see also Fig. 11 


The fluctuations of the occupations of the collective 
modes cause fluctuations in the qubit frequency and 
hence dephasing. For the qubit-cavity coupling, the so- 
called photon shot noise dephasing has been investigated 
in the 3D transmon architecture [34M35j . In particular, in 
Ref. IM] good agreement between theory and experiment 
was found over a range of average occupation number in 
the cavity h from small to relatively large (h ~ 3), and a 
residual occupation of order 1% was estimated. Since the 
collective modes are (weakly) coupled to the cavity and 
the residual occupation is small, here we restrict our¬ 
selves to the relevant case of small occupation number 
^ 1 for any collective mode. In this case each mode 
contributes a rate 1551 




> + 4x' 


with Kfi the decay rate of mode ^ 
dephasing rate 


= Er.. 


(70) 

], to the total qubit 

(71) 


Equation (70) was derived assuming that the effect of 
each mode can be treated independently [5S]. This ex¬ 
pression enables us to put an upper limit on the de¬ 
phasing rate, since independently of the rate satisfies 
bp < Ixplffp- The inequality is saturated for k,^ = 2|xp|, 
and becomes a strong upper bound in the limiting cases 
of /vp much bigger or smaller than |xp|. In the remain¬ 
ing of this section we restrict our attention to flux being 
zero or half a flux quantum, since it is experimentally 
established m that away from these “sweet spots” the 
fluxonium dephasing rate is determined by flux noise. 
However, we note that near resonances where the disper¬ 
sive shifts are enhanced, see Figs. and they could 
give rise to reproducible suppressions of coherence time 
T 2 at specific flux values. 

To determine if the collective modes can at least in 
principle be a significant source of dephasing, let us con¬ 
sider a worst-case scenario in which each rate Fp attains 
its maximum value and all collective modes are equally 
populated, = h. Then we have 




(72) 


We have calculated \xp,\ at zero and half flux quantum 
for both parameter sets in Table |lj see Fig. El Sum¬ 
ming over all modes and assuming n = 0.01, we arrive 
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FIG. 11. Absolute values of the total dispersive shifts |Xp| 
for parameter set 1 (inset) and 2 (main panel). Filled cir¬ 
cles are evaluated at zero external flux, empty circles at half 
flux quantum. The higher collective modes have very simi¬ 
lar values of |XpIi dominated by Xp^^ of Eq. (59). The large 
fluctuations of |Xp| for the low even modes are due to Xp of 
Eq. (39). In all cases, the lowest even mode (p = 2) has the 
largest shift. 


at the results summarized in Table El We find that in 
the worst case, the collective modes could limit the de¬ 
phasing time to about O.l^s at zero flux and Ifis at half 
flux quantum. Measured coherence times are longer than 
these estimates mm, indicating that the worst case is 
not realized in practice. In fact, in Ref. [TT]the Purcell- 
limited lifetime of the qubit at zero flux was measured to 
be at least 10/xs; then the results of Sec. IV C and Fig. 
indicate that the decay rate of the lowest even mode is 
of order 100 Hz, much smaller than the dispersive shift, 
and all other modes have even smaller decay rates IM- 
In the more realistic limit Kp ^ IX/jIj from Eqs. (70)-(71) 
we find 


r5 = h^«p, (73) 

p 


from which we get the estimates in the last column of 
Table corresponding to a dephasing time of order 1 s 
(dominated by the relaxation rate of the lowest even col¬ 
lective mode). This time scale is much longer than the 
coherence times measured in experiments, indicating that 
most likely the collective modes are not causing any sig¬ 
nificant dephasing. We caution the reader that the esti¬ 
mates in the last column of Table [H] rest mainly on the 
assumption (valid within our model for typical parameter 
values as in Table |T| that the collective modes are much 
more weakly coupled to the cavity than the qubit is, as in 
Fig. [3 if the assumption is not correct, this could result 
in a dephasing time shorter by several orders of magni¬ 
tudes - see also the end of Sec. Eml 
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r^(o) 

r^(o.5) 

T-'K 

i 

set 1 
set 2 

13.1 

7.99 

0.56 

1.17 

1 X 10"*’ 
1.5 X 10"® 


TABLE II. Estimates (in MHz) for the dephasing rates E^, 
Eq. (721, and T^, Eq. (731 - see text for details. 


VII. BROKEN PT-SYMMETRY: AN EXAMPLE 


In all the previous sections we have assumed the sys¬ 
tem to be PT-symmetric, which ensures the decoupling 
of qubit and odd collective modes in the (approximate) 
quadratic Lagrangian. In practice it is difficult to fabri¬ 
cate a perfectly symmetric circuit, so it is interesting to 
investigate what are the main qualitative consequences 
of breaking parity symmetry. To this end, we consider a 
simple case in which the symmetry is broken by taking 
the two coupling capacitors to have different values: 


= C, + CtSc/2 , C" = a - CtSc/2 (74) 


with Ct defined in Eq. (12). Therefore in this section Cc 


represent the average of the two coupling capacitors, and 
we have introduced the dimensionless asymmetry param¬ 
eter 


Hamiltonian H takes the same form as in Eq. (32): 


N-l 

P = ^ + Pi 


int 


Hv , 


P^ = - E*} cos ((/) -b <Pe) + ^El(, 


1 


ai:2 


H, = AEc,,p^^ + -EX, 

N-l 

Hint — ^ ^ dfiPfiPcj} ) 

(1=1 
N-l 

Hv = -'^ g^Pfj. eV 

/i=i 


(Jc)= 


4P^ 8Et 2Pg 


(80) 

(81) 

(82) 

(83) 

(84) 


-E^cP^eV 




Et 


E' 


Ec. 




2N 

M \-^9 


{E^f 4 


Despite the formal similarity, there are important dif¬ 
ferences between Eqs. (32) and (80): first, all collective 


modes appear in P, not just the even ones; in fact, we re¬ 
mind here that Pc./x is given by either Eq. (25) or Eq. (28) 


depending on the mode parity. Second, due to the asym¬ 
metry the qubit charging energy Eq is renormalized from 
the definition in Eq. (30): 


6c = 


gg - 

Ct 


(75) 


J_ _ J_ _ {6c)^ 
Pg ~ P^ 4Pt 


(85) 


with |Jc| < 1. Inclusion of the asymmetric capacitive 
coupling in the circuit Lagrangian amounts to the sub¬ 
stitutions [5H] 


Goo Goo ~ -TTv (<^c) , 

4Et 


Gofi —t G( 


0/T 


Sc 


-fi'-'ii 


2P“ t/msl 


(76) 

(77) 

(78) 


in Eqs. ([^-([^, anc 




8P 


Vmsl 


4eV (79) 


in Eq. (17). Note that matrix is not affected and, as 
discussea in Sec. [ml a rotation is needed to diagonalize 
it in the odd sector. This rotation in principle modifies 
the new term introduced in Eq. (©; we neglect this 
modification as it does not introduce any new qualitative 
feature - this is also a quantitatively good approximation 
if the parameter A, Eq. (13), is sufficiently small. 


Within the same approximations used previously [in 
particular, we assume again Eq. (31) to hold], the total 


Third, the coupling constants are different for even 
and odd modes: 


g^i = 


— NiNNlX. _ o^^c) . (86) 






The structure of Hamiltonian P in Eq. (80) shows that 


the main consequence of breaking the parity symmetry 
is the introduction of coupling between qubit and odd 
modes with coupling strength linear in the asymmetry 
parameter 5c. Therefore for strong asymmetry, |(5c| ^ 1, 
the odd modes influence the qubit in the same way as 
the even ones. Even for moderate asymmetry, [dcj ^ 0.1, 
the effect of the lower-energy odd modes may be non- 
negligible (at least near zero flux, where the qubit fre¬ 
quency is closer to those of the collective modes): while 
5c suppresses the coupling of the odd modes to the qubit, 
the odd modes with index 2p — 1 are closer in frequency 
to the qubit than the even modes with index 2p, and the 
smaller frequency difference generally increases the dis¬ 
persive shift, see Eq. (39); also, the term proportional 


to 1 /Eq in Eq. (84) roughly compensate for the 5c sup¬ 
pression of coupling between odd modes and cavity, thus 
giving similar lifetimes for odd (2p — 1) and even (2p) 
modes. On the other hand, small asymmetry at the per¬ 
cent level, as usually present in nominally symmetric de¬ 
vices, implies that the odd modes can be safely neglected. 
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FIG. 12. (Color online) Comparison between collective modes 
frequencies measured in Ref. [18] (red circles) with those cal¬ 
culated using the present theory (black crosses). The inset 
shows the ratio between theoretical and experimental frequen¬ 
cies. 


A. Comparison with experiment 


Asymmetrically coupled systems similar to that de¬ 
scribed above have been recently probed experimentally, 
which enable us to test in part our theory. For example, 
in Ref. [IHlan array of 80 junctions was placed in parallel 
to a (resonator) capacitor and the nine lowest resonant 
frequencies were measured . T he system is described by 
the Hamiltonian H in Eq. (80) if we set E^j = 0. Then H 
describes a harmonic oscillator linearly coupled to — 1 
oscillators. The resonant frequencies of the correspond¬ 
ing N independent oscillators can be easily calculated 
numerically; to compare with experiments, we note that 
the lowest mode in the experiment correspond to what 
we call the qubit mode (p, and therefore the higher even 
indices correspond to our odd modes and viceversa, odd 
mode indices in the experiments correspond to our even 
modes. 


In calculating the resonant frequencies we use as input 
parameters the array junction capacitance Cj = 40.1 fF 
{Eq ~ 483 MHz) and its Josephson inductance Lj = 
1.94 nH {Ej ~ 84.3 GHz), the capacitance to ground 
Cg = 114 aF {Eg ~ 170 GHz), the two coupling capac¬ 
itors C° = 5.5 fF and Cf = 3.16 fF {E^ ~ 2.24 GHz, 
dc ~ 0.27), and the resonator capacitor Cj = 3.19 fF 
{Eq ~ 6.07 GHz). The frequencies calculated with these 
parameters differ by less than 1% from the measured fre¬ 
quencies, see Fig. We note that while the array junc¬ 
tion parameters Cj and Lj agree with those reported in 
Ref. [iHl the ground capacitance Cg we estimate here is 
about three times bigger; since the calculations in the 
previous sections are based on the original estimate of 
Ref. [in they could underestimate, e.g., the dispersive 
shifts by almost one order of magnitude. 

A very recent experiment |28j reports the measurement 


of 14 resonant frequencies in an array of 200 junctions 
without shunting capacitor (Cj = 0). We can again 
compare our calculated frequencies with the measured 
ones: setting Cg = 98 aF and optimizing the other pa¬ 
rameters, we find again differences of less than 1% except 
for the third mode, whose measured frequency is about 
7% higher than the calculated one; this larger difference 
is likely due to the presence near the frequency of that 
mode of a spurious resonance |28j . 


VIII. COUPLING INTO THE 
SUPERINDUCTANCE 


So far, both for the PT-symmetric and the broken- 
symmetry cases, we have taken the coupling capacitors to 
be connected to the two islands separated by the phase- 
slip junction. The coupling capacitors can be attached to 
any island in the circuit, and in fact such a setup has been 
used in more recent experiments |11LI39] . In general, arbi¬ 
trary placement will immediately break parity symmetry 
even if the capacitances are the same for both capacitors. 
Here we consider briefly the simplest case of equal capac¬ 
itors placed symmetrically with respect to the phase-slip 
junction, so that parity symmetry is preserved and no 
qubit-odd mode interaction is allowed. Concretely, we 
take the first capacitor to be connected to island 6 > 0, 
where islands 0 and N are the two island surrounding 
the phase-slip junction; then the second capacitor is con¬ 
nected to island N — 6, and the maximum possible 6 is 
Sm = L(-^ ~ 1)/2J- In this configuration, the coupling 
Lagrangian Cy [cf. Eq. 0] becomes 


Cv 




(peV + 



S2p 


VpeV. 


(87) 

This formula correctly reduces to Eq. ( |T7| for (5 = 0, 
while for 5 > 0 a new coupling between cavity and even 
modes is present. 

Changing the position of the coupling capacitors also 
affects the kinetic energy part To of the Lagrangian [cf. 
Eq. ([^], and a general treatment of this modified term is 
quite cumbersome. Here we consider the simple limit in 
which we neglect the ground capacitances, Cg,Cg —0. 
In this case, as we show in Appendix]^ the Hamiltonian 
is 


H = H^ + H2 + Hint +Hv + Y,Hg, (88) 

~ > ( 89 ) 

^2 = 4 + lEje2 , ( 90 ) 

gJ^EQ2 , ^ 

^ -13^2-P2P0 = g2P2P<l> , (91) 

Hg = AE^qpI + ]^E<}ep , 


(92) 
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FIG. 13. (Color online) In all four panels, filled circles are 
used to denote parameter set 1 and empty ones for set 2. 
The horizontal bottom (top) scales are used for set 1 (2). a) 
Normalized effective qubit charging energy Ef^{5)/Efg(Q) vs. 
coupling capacitor position 5. b) Normalized collective mode 
charging energy i?2,efi((5)/i52,eff(0) vs. 5. c) dimensionless pa¬ 
rameter g[Eq. (96l] vs. 5. d) qubit-collective mode coupling 
g 2 [Eq. (1^] vs. 5. 


Hv = - 


+ 



Ap2eV Ec2 
l-g2^ 

= -g4>,cP(j,eV -I- g2,cP2eV 


(93) 


where the parameters g 2 , 
well as 


3 , and g 2 ,c defined above as 
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Ec2 E‘. 
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2E: 
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26 

N 
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25 

iV 


g' 




2{E^f V N 


(94) 

(95) 

(96) 


depend on the position 5 of the coupling capacitors. 

The main qualitative feature of the Hamiltonian in 
Eq. (88) is that the qubit and the cavity both couple to 


only one collective mode whose charging energy is renor¬ 
malized below E^, while all the other N—2 modes remain 
degenerate and uncoupled. (Of course in the presence 
of ground capacitances the degeneracy is lifted and all 
the even modes couple to both qubit and cavity.) While 
for typical experimental parameters the effective qubit 
charging energy E^^ = Eq/{1 — g^) moderately increases 
as 5 increases towards N/2, the collective mode effec¬ 
tive charging energy if 2 ,eff = Ec 2/(1 — g^) can be more 
strongly suppressed when 5 ~ N/4, see Figs. 13 1 and 13 d. 
The dimensionless parameter g varies non-monotonically 
as function of 5 and is generally small, see Fig. [T^. The 
qubit-collective mode coupling strength g 2 also depends 
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FIG. 14. (Color online) Normalized qubit-cavity coupling 
[Eq. (93l] vs. coupling capacitors position 5. Horizontal bot¬ 
tom (top) scale is used for set 1, filled circles (2, empty circles). 
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FIG. 15. (Color online) Collective mode-cavity coupling g 2 ,c 
[Eq. (93l] vs. coupling capacitors position 5. Horizontal bot¬ 
tom (top) scale is used for set 1, filled circles (2, emp ty c ircles). 
The squares at 5 = 0 are given by gi/iE^ [cf. Eq. (36l], with 
gi from Eq. (351, and are plotted here for comparison. 


significantly on 5, see Fig.[T3p; note that the largest val¬ 
ues of g 2 at 5 ~ N/8 are a significant fraction of (or 
comparable to) the coupling strengths between qubit and 
lowest collective modes calculated for 5 = 0 but in the 
presence of ground capacitances, see Fig. These ob¬ 
servations imply that by appropriately placing additional 
capacitors in the array, both the collective modes spec¬ 
trum and the coupling strength with the qubit can be 
controlled to some degree. In particular, to minimize the 
effects of the collective modes on the qubit the coupling 
capacitors should either be placed next to the phase-slip 
junction (5 = 0), or opposite to it (5 = 5m), while inter¬ 
mediate positions (especially in the range 5 ^ N/S—N/A) 
maximize those effects. 

In Eq. (931 we give expressions for the qubit-cavity and 
collective mode-cavity dimensionless couplings g^^c and 
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g 2 ,c- For (5 = 0, reduces to the value in Eq. (36); as 5 
increases, it gradually decreases down to a value approx¬ 
imately 1/iV times the initial one when 5 = 6m- In con¬ 
trast, g 2 ,c takes its smallest values for (5=1 and (5 = Sm, 
where it is approximately given by g™” « AEq/^/2Eq; 
the largest values at 6 ~ iV/4 are bigger by a factor of 
less than 3. It turns out [cf. Fig. 15 that for typical ex¬ 
perimental parameters even the minimum value g™“ is 
larger than the strongest collective mode-cavity coupling 
gi/AE^ in Eq. (36) |40) . 

The contrasting dependence on 6 of the two couplings 
indicates that moving the coupling capacitors away from 
the phase-slip junction can adversely affect the qubit co¬ 
herence: indeed as 6 increases, the coupling capacitors 
must also increase to attain the desired qubit-cavity cou¬ 
pling strength; both moving the capacitors and increas¬ 
ing their capacitance, however, raise the collective mode- 
cavity coupling, which in turns increases the collective 
mode (Purcell) decay rate k (cf. Sec. IV C) and hence the 
qubit dephasing rate, see Sec. lYll We therefore conclude 
that placing the coupling capacitors beside the phase-slip 
junction, (5 = 0, is the optimal choice. 

If the coupling capacitors are placed opposite to the 
phase-slip junction ((5 ^ N/2), the coupling capacitance 
should be increased by a factor of order N to compensate 
for the decrease in qubit-cavity coupling strength. To¬ 
gether with the stronger mode-cavity coupling, g™ as 
compared to gi/iE^, this increase would raise the col¬ 
lective modes decay rate by at least ^ lOiV^ (i.e., ~ 10"^). 
Then the decay rate of the lowest even mode would be 
faster than that of the qubit, the dephasing rate F^ of 
Eq. ( [T^ would also increase by about 4 orders of mag¬ 
nitudes, and the corresponding dephasing time would be 
about 0.1 ms. This time is one order of magnitude longer 
than the coherence time measured in Ref. 1111 indicating 
that the collective modes are not limiting coherence in 
current experiments. On the other hand, our estimate 
is much shorter than the measured ~ 10 ms relaxation 
time at half flux quantum, so the effect of the collective 
mode could in principle be observable, if other dephasing 
mechanisms can be identified and suppressed. 


IX. SUMMARY 


Fig.[g 

In Sec. we consider the leading-order non-linearity 
of the array junctions, which introduces additional qubit- 
collective mode interactions. Among these interactions, 
the term which leads to the strongest dispersive shifts 
effectively induces fluctuations in the qubit inductive en¬ 
ergy when the collective modes are excited, see Sec. |VB 
As we discuss in Sec. |VI[ the total dispersive shifts (i.e., 
including also the effect of ground capacitances) are much 
bigger than the collective mode decay rates, so the lat¬ 
ter determine the qubit dephasing rate. We find that 
the collective modes do not significantly contribute to 
dephasing, so long as they are more weakly coupled to 
the cavity than the qubit is; the weak coupling is generi- 
cally achieved if the qubit-cavity coupling capacitors are 
placed next to the phase-slip junction. However, we es¬ 
timate in Sec. IVIlIl that the collective-mode induced de¬ 
phasing could become observable if the coupling capaci¬ 
tors are placed opposite to the phase-slip junction. 
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Appendix A: Lagrangian 


This appendix provides a derivation of the Lagrangian 
£ for the collective and qubit modes [Eq. (18)] starting 
from a more familiar “textbook” formulation in terms of 
phases and voltages. To this end, we split £ as a sum of 
kinetic and potential energy parts as usual, C — T — Us, 
and write the potential energy as 


N 

Us = —Ej ^ cos 9m — Ej cos j ^ 9„ 


N 


£e 


(Al) 


m—l 


\m—l 


In this paper we study the collective modes in the ar¬ 
ray of Josephson junction forming the superinductance of 
the fluxonium qubit. We derive an approximate Hamil¬ 
tonian, Eq. (32), that includes the interactions between 
the qubit mode and the collective modes in the presence 
of ground capacitances. The approximations place some 
restriction on the number of array junction to which the 
model applies, see Eq. (311, but this condition is in prac¬ 
tice much weaker than that given by the array “screening 
length” [Eq. ( [^ ] and it is satisfied in current experi¬ 
ments. A generalization of this Hamiltonian enable us to 
favorably compare the calculated spectrum of the collec¬ 
tive modes to two recent experiments, see Sec. VH and 


Here 6m is the (gauge-invariant) phase difference across 
junction m in the array, and the first term on the right 
hand side is the Josephson energy of the array junctions. 
The last term in the above equation is the phase-slip 
junction energy, and in writing this term we have taken 
into account the fluxoid quantization condition 

N 

'^0m+^e = 2?™ (A2) 

ttl—Q 


with n integer. 

The kinetic energy part T is more easily expressed in 
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terms of the voltages of each island: 


T = Ts + Tg + Ty , 


m—1 

N 


i6£;g, 


16i?^ 






i6£;r ’ 

m—0 9 


rp {‘Pm — 2eVm) 

Tv= 


771=0 


16^;? 


(A3) 

(A4) 

(A5) 

(A6) 


Here Tg is the charging energy due to the junctions 
capacitances, To due to capacitances between each is¬ 
land and ground, and Ty due to capacitive coupling to 
external voltage sources Vm- In the above equations 
= e^/2C™ and = 6^72(7™ are charging ener¬ 
gies of ground and coupling capacitors for the mth island 
and they can in general be different for each island. We 
stress that all equations in this Appendix are valid for 
this generic case, not just for the specific circuit depicted 
in Fig.[lJ 

To rewrite T in terms of the phase differences 0m, we 
use the relationship 


Pm — PQ + 1 (A7) 

1=1 


valid for m = where we have taken as a 

reference phase. Then it is straightforward to write Tg 
in terms of 6m variables 


N hi 

Tg = V 
* ^ \QE% 

771=1 ^ 


(S 


N 

771=1 * 




(AS) 


The other two terms in T take the form 


•2 ^ 1 
tg=^.ey: 


Tv = 


16 £;o ' 16 A" 

9 771=1 9 

{po - 2eVof 
16E0 
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^0 


z=i 


(A9) 

(AlO) 


+ E + XI “ 2eI4, 


771=1 


Z = 1 


manipulations lead to 
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1 ^ 

Tg = X SmnSmOn , 


(All) 


771,71=1 

min{r7i,7i} —1 ]\[ 


5nm-Et X \Ei ^ E’-ji ^ E^r 

i=0 i=max{m,n} ^ 9 J 
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771=1 7=771 
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where 
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i=0 


El 
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(A13) 


These equations correctly reduce to those of Ref. HH in 
the absence of coupling capacitors. 

As a final step, we introduce a new set of variables via 
the relations 


N 


4’ — 'y ^m , 

771=1 

N 

= E 


(A14) 

(A15) 


with index p, = N — 1 and inverse 6m = <7’/-A -I- 

JYu. W^m4.fj.- The matrix W^m must satisfy the conditions 


W^mW,, 


= s„ 


In terms of 


Em = 0 and Em , 
these new variables we find Tg as in Eq. (0 and Ug as in 
Eq. ^ . Formulas for Tq and Ty and arbitrary kF^m are 
not instructive, so we do not report them here. For the 
specific circuit configuration and choice of IF^m described 
in Sec. |Tl| the corresponding formulas are given there 
and follow directly from the equations above. Modifica¬ 
tions of those formulas for a different circuit configuration 
breaking parity symmetry are discussed in Sec. Ivnl A 
third circuit with coupling capacitors connected into the 


array is briefly considered in Sec. VIII Here we mention 
a useful identity valid for the choice of Wf_im in Eq. ([^ : 


E 


1 


(A16) 


[for the notation used, see the definitions in Eqs. (14)- 

(del)]- 


We next note that C is independent of po, so that 
d£/dpo is a conserved quantity, the total charge of the 
circuit [22]. Using this conservation law we can express 
Po in terms of the variables 6m and thus eliminate it 
from the Lagrangian. In this way, standard algebraic 


Appendix B: Choice of parameters 


An often measured property of a flux-tunable qubit 
such as the fluxonium is its spectrum as a function of 
flux, u}io{f). The spectrum can be obtained by numerical 
diagonalization of the qubit Hamiltonian E[^, Eq. (33), 
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where the inductive energy should be replaced by 
[Eq. ([ 5 ^]. For the experiments reported in Refs. Eland 
mi with = 43 and TV = 95 array junctions respec¬ 
tively, this procedure leads to the parameters reported 
in Table |III| We also give there our rough estimate of 
the ratio Cg/Cg which is based on the geometry of the 
devices in the two experiments. 


The phase-slip junction Josephson energy Ej in Ta¬ 
ble |l] is taken directly from the experimental estimates 
in Table [TT^ For the coupling capacitor energy Eq, we 
use for set 1 the value of coupling capacitance given in 
Ref. El while for set 2 we take as an example the value re¬ 
ported in Ref. 1181 The latter experiment was performed 
with junctions fabricated with the same procedure used 
to fabricate the fluxonium junctions of Ref. [11] for this 
reason, we also use the value of ground capacitance in 


Ref. 


I to estimate Eg for set 2. The value of Eg for 


set 1 is smaller than that of set 2 by a factor of 0.4 be¬ 
cause, according to the supplementary to Ref. UHl the 
different fabrication processes lead to such a difference 
in capacitances to ground. Values for Fg then follows 
from the geometrically estimated ratio Cg/C®. To fur¬ 
ther constraint the parameters, since the only difference 
between phase slip and array junctions is in their area, 
we assume that their plasma frequency is the same, so 
that EjEq = E^jEq. We then choose t he p arameters 
Ej, Eq, and Eq as to obtain, using Eqs. ( [^ and (53l, 
the values of Eq and Eq reported in Table]lll[ 


To eliminate the off-diagonal terms at order A, we take 

51 to have elements Si^pp = 0 and for p 7 ^ cr 

Si,p. = ^ , (C2) 

-^d,<7(7 -^d,pp 

where for example the elements TajPo- of matrix T\ are de¬ 
fined via Tx = J2p o-^a.pctCpCct- We can similarly choose 

5 2 to eliminate the off-diagonal terms at order A^. The 
remaining diagonal elements are then, up to order A^, 


Td,pp + XTx,pp-^[Si,Tx]gg . 

The last term is more explicitly written as 

TxnaTx,ap 


ia.ui„ = 2i: 


a^p 


Td,a<j — Td, 


(C3) 


(C4) 


pp 


By substituting the matrix elements into Eq. (C3), we 
arrive at Eq. (251. 


Note that the rotated modes C,p are related to the orig- 
inal modes Cct via Cp = (s'®) ^C<t- For use in Sec. 


VA 


we define the matrix that performs the rotation m 
the odd sector while leaving the even sector unchanged: 


Ap^ 


p, i/even (C5) 

0 otherwise 


Appendix C: Rotation in the odd sector 


The diagonalization of the kinetic energy matrix To 
[Eq. ( |^ ] of the odd sector can be seen as a rotation 
among the odd-sector coordinates Cp- Such a rotation 
can be constructed perturbatively order by order in the 
parameter A, with a procedure analogous to that used to 
perform a Schrieffer-Wolff transformation to an effective 
Hamiltonian: we want to obtain an antisymmetric matrix 
S = J2n=i such that the the product e~^Toe^ is 

diagonal. Expanding this product up to second order we 
have 

e-^Toe^ =Td + XiTx-[Si,Td]) 

+ A2 Q [5i, [5i,Td]] - [S2,Td] - [Fi,rA]) + ... 


Appendix D: Even sector Hamiltonian 


The Lagrangian Ce in Eq. ( [27| has a quadratic kinetic 
energy part, so that transforming to the Hamiltonian 
amounts to inverting the {No -|- 1) x {No + 1) matrix 
K = Kd + Ki whose diagonal part Kd has entries 


Aid,00 — 

^d,pp — 


8e: 


c 


8E: 


c,p 


p=l,...,No 


and the interaction matrix Ki has elements 


(Dl) 

(D2) 


— ^i,pQ — 


1 C2p 


32VmE-slg 


(D3) 


and all other elements are zero. Indeed, after performing 
a standard Legendre transformation the kinetic part Tlx 
of the Hamiltonian has the form 



N 

A) 

Efi 

El 


set 1 

43 

8.93 

2.39 

0.52 

32 

set 2 

95 

10.2 

3.60 

0.46 

0.6 


TABLE III. Additional fluxonium parameters estimated as 
explained in AppendixEnergies are given in GHz. 


Hk = — nK , (D4) 

where 

n = (n0,ni,... ,n7vj , (D5) 

u=(-eV/4F^,0,...,0) (D6) 
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are two A^e + 1-dimensional vectors. We can formally 
write the inverse matrix as 


K-^=K7^ 


1 


1 - K,K-^K,K-^ 


— K~^K K~^ 

j\d 


1 - K,K-^K,K-^ 


and approximate it as 
K- 


Kr-Ka^K,K- 


(D 8 ) 


provided that entries of matrix KiK^^KiK^^ are small 
compared to unity. Given the structure of matrix Ki, 
this condition translates into 




^C,p^C C2p C2a 

32N {E^f4p sL 


< 1. (D9) 


If the perturbative condition Eq. ( [^ is satisfied, we can 
approximate Eq ^ « Eg, - see Eq. (28); using that in 


order of magnitude E^ ~ Eg, it then follows that the 


condition (D9) is also satisfied. More interesting is the 


opposite regime in which Eq. (23) is violated; then we can 
take E^ p « 4E“s|^. Substituting this approximation 
into Eq. (D9) and expanding the trigonometric functions 


for small p and a (which maximizes the left hand side), 
we arrive at Eq. ([M|) . Using Eq. (|D4|) and (|D 8 |) we arrive 


at Eqs. (32|-(36). 


In this way we find 

.Ueff — Hfj) 


(D7) with 


1,3 


- ^ [E[p + El^p\ 


~ gp 

4 


(E5) 
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e,- - ei 


■ 2uju 


nn -b 


(E 6 ) 


When projected onto the qubit subspace, the last line 
above gives a small contribution to qubit frequency shift 
as well as to the dispersive shifts; note that due to our as¬ 
sumption wio < u!p, the terms in the last line are always 
finite. To see their smallness, consider the inequalities 


E 


9p0j 
G ~ Co + 2a;n 


± ' 


glu 


e, - ei 


2u!ii 


< 


E 


9p-0j 


glu 


ELil 

8N 


OJp 

4 TP, c"2 




W, 


2 


(E7) 


In the last step we used that since state |0) (|I)) is mostly 
localized at potential wells with minima between 0 and 
~ ±7r (~ ±27r and ~ ±7r) we have 


Appendix E: Schrieffer-Wolff transformation for u: 


( 2 ) 




In Sec. VB we have argued that the the part propor¬ 


tional to of the potential term gives rise to 

fluctuations in the inductive energy and hence to corre¬ 
sponding dispersive shifts; here we consider the remain¬ 
ing part of proportional to -I- apUp. 

^^ = EE gpij\l){j\ (ajjaj) -b a^ap) , (EI) 

M l,j 


where 


gpij — 


E‘]ll ( 6\ ElII 

^(^1 1-cos^ b-)-^(^|</>^b'). (E 2 ) 


An effective qubit-collective modes Hamiltonian Ees 
that includes the effect of this interaction term can be 
obtained by keeping the next-to-leading order, diagonal 
part of e~^He^, where 


H = H^ + J2Hp + U2 


(E3) 


ianoi<[a + iM 


Z = 0,1. 


(E 8 ) 


The last expression in Eq. (E7) shows that the dispersive 


shift from the last line of Eq. (E 6 ) is much smaller than 

see Eqs. (59) and (61). _ 

The term on the first line in Eq. (E 6 ) can diverge if 


the resonance condition ej — ei — 2ujp = 0 is met, but 
this divergence simply signals the breakdown of the per¬ 
turbative approach when — £/ — 2u>p « g^ij- There¬ 
fore the coefficient of the largest possible contribution 
from the first line of Eq. (E 6 ) is much smaller in magni¬ 


tude than Tr^ELip/2N (the smallness is due to the matrix 
element begin evaluated between one of the low-energy 
qubit states, | 0 ) or | 1 ), and a state with much higher en¬ 
ergy). Hence we find that even near resonance this term 
is smaller than ) moreover, one should keep in mind 
that the latter shifts the qubit frequency if any collec¬ 
tive mode has at least one excitation, while the former 
changes the qubit frequency only if a specific mode (the 
one near resonance with a qubit transition) is excited at 
least twice, and the probability for the latter situation to 
happen is smaller by a factor ~ fi/N if n is the average 
occupation probability - see also the next Appendix. 


with of Eq. (33) and Hp of Eq. (82), and 


S = EE^^blOO’l 

M 1,3 






e, - £/ - 2ujp e, - ei + 2u!p 


(E4) 


Appendix F: Occupation probabilities 

In this Appendix we comment briefly on the assumed 
smallness of the occupation probabilities for the collective 
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modes. We assume for simplicity an equilibrium proba¬ 
bility for each mode, so that the probability of having 
Up excitations in mode p is 


= 


1 


1 -I- rip \ 1 -I- rip 


(FI) 


where hp is the mode average occupation. For an order 
of magnitude estimate, we take the latter to be the same 
for all modes, hp = h, and small, h <C 1. Then the 
probability Pq that none of the — 1 collective modes is 
occupied is simply 

^0= ^l-(iV-l)h (F2) 

[l + n) 

while the probabilities Pi that one mode has a single 
excitation, P 2 that one mode has two excitations, and 
Pi 1 that two modes have one excitation are 


Pi = (lV-l)- -(F3) 

(1 + h)^ 

_ 2 

P2 = (iV-l) (F4) 

{l + ny ^ 


(l + n) 


(N-l)^ 


n\ 


In approximating the above formulas, we assumed h <C 
1 /A^ ^ 1, and it is evident that under these assumptions 
we have P 2 <C Pip < Pi < Pq. 


Appendix G: Derivation of potentials and 
and calcnlation of freqnency change Aiaio,^ 

The starting point to derive the formulas for and 


/o\ 

[Eqs. (62) and (67), respectively] is the following 
approximation 
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sin 


Wprn^p 




6 


m,/i, 1^,(7 


(Gl) 

where we used the property Wpm = 0 to elim¬ 
inate the lowest order contribution. Next, we ex¬ 
press the collective mode coordinates in terms of cre¬ 
ation/annihilation operators and after normal ordering 
we find 






2V2 


where H.c. denotes the Hermitian conjugate. The last 
term in square brackets gives the linear interaction term 


■ 


1 


^ E E Cm ■ (G3) 

m 

To proceed further, we note that 

y] = I/Vm (G4) 


if one of three conditions is satisfied: 

p + u — a = 0, 
/r — + cr = 0, 

p — v — a = 0, 

while 


(G5) 

(G6) 

(G7) 


y] WpmW^mW^m = -1/V2N (G8) 


if 


p + 1/ + (T = 2N . 


(G9) 


The sum over m vanishes otherwise. For the sum in 


Eq. (G3) this only leaves two possibilities, p = 2v and 
p = ^ — v) , which imply that p must be even. Setting 
p = 2p and using the definitions o f rjp and Ng given after 
Eq. (27), we can finally write Eq. (G3) in the form given 
in Eq. (62). 


Using again the identities in Eqs. (G4)-(G9), it is 


straightforward to cast the terms with products of three 
operators in Eq. (G2) in the form of Eq. (67). Here we fo¬ 


cus on the calculation of Awio,ii, the change in frequency 


( 3 ) A Iv 

due to the interactions in , Eq. (67), when a single 
collective mode is excited. To this endTlet us define the 
corrections (5e/,o and (5epp to the qubit energy e; depend¬ 
ing on whether the collective modes are in their ground 
state or in state jlp), respectively. These corrections can 
be calculated at second order in perturbation theory in 


for example: 




Se. 


1,0 


= E 

jj 


e; - G - Ef 


(GIO) 


Here jj, /) denotes a generic state jj) for the qubit, with 
energy ej, and some Fock state ]/) with energy Ef for 
the collective modes. The similar expression for epn is 




3^1,p = E 


jJ 


ei+ujp- ej - Ef 


(Gll) 


In terms of these corrections the qubit frequency change 


[{alalal + Sa^alaa -k H.c.) + 3 (aj, -k Qp) S^a] , 

(G2) 


^wio,p — — Seo,p) — (i5ei_o ~ <5eo,o) • (G12) 
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A great simplification in calculating Awio,;i is achieved by 
noticing that in the difference dei^p,—dei g all contributions 
originating from terms in Eq- (67), for which none 
of the indices of the operators in thar^uation coincides 
with fi cancel out. In other words, only terms for which 
at least one index is p. can contribute to the frequency 
change. Moreover, since the Hermitian conjugate terms 
in Eq. (67) contain two or more annihilation operators, 
they also do not contribute to Awio,^. To conc retely cal¬ 
culate the matrix elements entering Eq. (Gil), we need 


to consider the action of the terms in square brackets in 


Eq. (67) onto the singly excited state |1^); for example. 


the first one gives: 




|l/jl/il;yl/i+y) (G13) 

+2V^lSplp+i/ljy) + y/2\2filp-i^l^) 

V V 

+V6|3^12m) + (v^)'|2^2p/2), 


where the prime at the summation symbols implies that 
all the collective mode indices in the states being summed 
must be different. The first term on the right hand side 
is the one in which no index coincides with p and, as dis¬ 
cussed above, this term does not contribute to the differ¬ 
ence 6ei^p, — Seifi. The last two terms arise from particular 
combination of the indices = v = p and fj, = v = p/2, 
respectively); the last one is present only if p is even. 
Similar terms in which two indices are equal also arise 


from the other operators in square brackets in Eq. (67) 


m- However, we discard these terms in comparison with 
the terms with sum over index v (see the second line in 
Eq. ( |G13 )): because of the sums, the discarded terms are 
smaller by a factor of order 1/A. Taking these consider¬ 
ations into accounts, standard calculation of the matrix 
elements for the collective modes gives 
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(G14) 


^2 p2 


- 1 - 

C/ Cj UJn UJ2N—Ti—i' Cl Cj -\- UJn UJ,^ UI2N — Ti—h 


where the approximate equality indicates that we are ne¬ 
glecting 1/A corrections originating from terms like the 
last two in Eq. (G13). 

To estimate a bound on the frequency change Awio.p, 
we note that for Z = 0 , 1 so long as Wio < 2 a ;i — u in-i 
there are no divergent contributions in Eq. (G14) and 
the second terms in round brackets are larger than the 
first ones. For an order-of-magnitude estimate, we fur¬ 
ther neglect the dependence of the oscillator lengths 
and collective mode frequencies on th e array ground 
capacitance and substitute in Eq. ( |G14 ) the values £o 
and Wp = yjSE/jE/j, respectively. In this way we find 
(for Z = 0 , 1 ) 
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Wp + (Cj - Cl) 


With the performed approximations, there is no depen¬ 
dence on collective mode indices, so that the sums in 
round brackets can be bounded by A. An upper bound 


for the term in square brackets is then 2 A/cli“ [32]. Fi¬ 
nally we approximate the qubit matrix element as fol¬ 
lows: 
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(jjAsin—|Z) 


\mYi)\ 


(G16) 


Rough upper bounds for the latter matri x ele ments are 
(27r)^ for Z = 1 and tt^ for Z = 0, see Eq. (E 8 ). We thus 


find that the bound is tighter for the Z = 0 correction 
compared to the Z = 1 one, and hence 

|Aa;io.p| < 2 \Sci,-, - 5q,o| ^ (G17) 


We now compare the frequency change Aa;io,p with 


that due to the dispersive shift x/i ^ originating from 


( 2 ) A k 

the quadratic interaction , see Eq. (59); the latter 
frequency change is given, in order of magnitude, by the 
upper bound in Eq. (61), so that an approximate bound 


on the ratio between the two quantities is 


|Aa;io^p| 

2 xf" 


< (^4)" 


El 
wio ’ 
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where we used E'j/ujp = £q^. For typical experimen¬ 
tal parameters, the right hand side is at most of order 
unity for any value of flux. Given that our approxima¬ 
tions place a loose bound on Awio.p, we conclude that 
the latter can generally be neglected in comparison with 

A/i 


This transformation is a rotation for any vector v with 
the definition 


= E^ 

(T — l 


(H 6 ) 


Appendix H: Derivation of Eq. (88) 


For the circuit considered in Sec. VIII (symmetrically 


placed coupling capacitors and no ground capacitance) 
the matrix G entering Eq. (1^ takes the form: 


Goo — 
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2E: 


Gofj. — 


2S 

N 


E: 
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1 - 






26 
fV 

1 ^2fi6 ^2iyS 

N Sa Si, 


S2^S 






0^J,+ 1 ■ 


(HI) 

(H2) 

(H3) 


Note that both matrix G and Cy in Eq. (871 vanish for 


5 = N/2] this choice (possible only for even N) connects 
both capacitors to the same island, and the vanishing is 
due to the arbitrariness in choosing a reference for the 
electric potential. For generic Q < 5 < N/2 we see that, 
beside a renormalization of the qubit charging energy 
[Eq. (Hll], the co upli ng capacitors couple the qubit and 
even modes [Eq. similarly to the effect of ground 

capacitances G“ in the array, while at variance with the 
effect of the ground capacitances, the coupling capaci¬ 
tors lead to mode-mode interaction in the even sector 


[Eq. (H3)], rather than in the odd one. (Of course par¬ 
ity symmetry still guarantees the qubit-odd modes de¬ 
coupling.) The situation, however, is only apparently 
complicated, since a change of variables shows that the 
problem reduces to a single collective mode coupled to 
the qubit and the cavity, while all other modes remain 
degenerate and decoupled. Indeed, let us introduce the 
new even-sector variables: 


N, 


E Cp 
p—1 ' ' 

klp-l ^ 


VpVg 

V\p\v\p-1 


-Va ■ 


(H4) 

(H5) 


for the norms, and |z)| = [i^lAfe- In the present case we 
take the components Vp to be 


1 S^pS 


^ ^/N S2p 

and find that the kinetic energy term Tq simplifies to 


(H7) 


16 


Goo</>" + - 2 


1 - 26/N 

V2E^ 


’\m<P 


(H 8 ) 


Similarly, £y of Eq. (87) becomes 


Gv = -tE (^-%)'PeV+ \v\^ieV . (H9) 


4:E^ V N 


The norm of vector v can be calculated explicitly 


J 1 - 


2(5 

iV 


(HIO) 


Neglecting the array non-linearities, the total La- 
grangia n is t he s um of Cu in Eq. ([^ with Tq and £y 
of Eqs. (H8)-(H9). Then performing the Legendre trans¬ 


form and defining 


c = J Vp /2 , M even 
' Cip+i )/2 , M odd 


(Hll) 


for /i = l,...,iV — 1, we arrive at the Hamiltonian in 
Eq. (ISSl). 
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